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Abstract. Varying Chaplygin gas is one of the dark fluids actively studied in modern cos¬ 
mology. It does belong to the group of the fluids which has an explicitly given EoS. From the 
other hand phase space does contain all possible states of the system. Therefore, phase space 
analysis of the cosmological models does allow to understand qualitative behavior and esti¬ 
mate required characteristics of the models. Phase space analysis is a convenient approach to 
study a cosmological model, because we do not need to solve a system of differential equations 
for a given initial conditions, instead, we need to deal with appropriate algebraic equations. 
The goal of this paper is to find late time attractors for the cosmological models, where a 
varying Chaplygin gas is one of the components of the large sale universe. We will pay our 
attention to some non linear interacting models. 
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1 Introduction 

In modern cosmology there are several key open problems and varions approaches to solve 
them [1] - [8]. In this work we have a goal to consider varions cosmological models and obtain 
appropriate conditions to have late time accelerated expansion of the nniverse [9] - [16]. Phase 
space portrait describing the dynamics of the model is one of the mechanisms allowing ns to 
have a qnalitative nnderstanding of the sitnation. This does allow ns to have all states of 
the model and we do not need to solve a system of differential eqnations, instead, we need 
to work with appropriate algebraic eqnations [17] - [25]. In this section we will provide some 
aspects related to the cosmological models which we will stndy in this work. Dark fluids are 
actively nsed to explain an accelerated expansion of the large scale nniverse. One of the dark 
flnids it is Chaplygin gas [26] - [29] 


Pc — Ape , (1-1) 

Pc 

where A^ B and a are positive constants and pc it is the energy density of the gas. In Physical 
Literatnre there are varions modifications of this flnid [30] - [36]. One of the modifications 
it is varying Chaplygin gas, where it is assnmed that either A oi B are not constants. It is 
possible that both of them are not constants. In this paper we will consider two models of 
varying Chaplygin gas given (see for instance [30] and [36]) as 


Pc — Ape 


BH-^ 
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Pc 


( 1 . 2 ) 













and 


Pc = Ap, 


Ba-^ 


(1.3) 


C 


Pc 


where H it is the Hubble parameter, a it is the scale factor, while n it is a constant. A 
phenomenological assumption about the content of the universe it is a usual procedure in 
modern cosmology. It is well known that to describe the dynamics of the large scale flat FRW 
universe (supported from observational data) we need held equations of General Relativity 


^2 ^ ^ ^ SttGp 


(1.4) 


a AttG , , 

- =-—(/J + 3P), (1.5) 

a 3 

where p it is the energy density of the effective huid, while P it is the pressure. Furthermore, 
we will work with the models, where the cosmological constant A = 0, the gravitational 
constant G and c are constants with c = SvrG = 1. It is easy to see from Eq (1.4) and Eq. (1.5), 
that we should make an assumption, for example, about the structure of the effective huid 
giving appropriate expanding behavior of the universe. From observational data it is known, 
that two kind of energy are involved into the dynamics of the large scale universe. This does 
mean that the energy density and pressure of the effective huid could be approximated as 


P = Pdm + PDE-, 


( 1 . 6 ) 


and 

P = Pdm + Pde, (1-7) 

where pde and Pde are the energy density and pressure of dark energy, while pdm and Pdm 
are the energy density and pressure of dark matter. There are various assumptions about 
the dark energy and dark matter, but all of them are still assumptions. There is another 
interesting fact, that should be mentioned about Chaplygin gas Eq. (1.1). It is a joint model 
of dark energy and dark matter. Some critics on this issue is also does exist in Physical 
Literature (see for instance [29]). There is a huge number of viscous dark huids [37] (and 
references therein). Dark huids could be either linear or non linear. Dark huids could be even 
solutions of different algebraic or differential equations [38]. Explanation of the accelerated 
expansion of the large scale universe could involve additional ideas despite to dark energy. 
In Physical Literature, there are various ideas on this topic. From the other hand it is well 
known, that modern cosmology has several other problems, and we believe that the solutions 
some of them are hidden into the other ones. An interesting situation has been noticed 
recently. Considered various cosmological models have proved that an interaction between 
dark energy and dark matter is a phenomenologically good idea. Practically there is nothing 
against to this idea. An interaction between dark components of the universe, can improve 
some unpleasant aspects of the models. This is a right motivation to continue a research 
on interacting cosmological models. Particularly, when there is an increase of the number 
of the interacting viable models. There are different parameterizations of the interaction 
term. Some linear and non linear interactions are used intensively. Recently, sign changeable 
interactions started to be used very intensively [39] - [62]. Interaction in this context should 
be understood as either dark energy transfers into dark matter or dark matter transfers into 
dark energy. If we consider interacting models, we do mean that the following should take a 
place 

Pdm + ^H{pdm + Pdm) = Q, 
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( 1 . 8 ) 





and 


PDF + ^H{pde + Pde) =—Q- (1-9) 

To construct cosmological models of this paper, we assume that dark matter is a pressureless 
fluid. The other aspect allowing us to study our models it is an assumption of the form of 
the interaction term Q. Particularly, we are interested in non linear interactions, which will 
be presented in section 2. 

The paper is organised as follows: In section 2 we will present basics on the phase space 
analysis to find late time attractor solutions of the field equations of General Relativity. 
In section 3 phase space analysis is performed, late time attractors are found and classi¬ 
fied according to their cosmological applicability. Finally, discussion on obtained results are 
summarised in section 4. 

2 Interacting models and autonomous system 

In Physical Literature there is a huge number of works devoted to phase space analysis of 
various cosmological models. In order to start phase space study of a model, an appropriate 
autonomous system should be found first. The critical points are solutions of the autonomous 
system. They are stable if appropriate Jacobian matrix has negative trace and positive 
determinant. This is in case of linear stability. Following Ref. [60] for our models we set 


Pc 

"" 3^2’ 

(2.1) 

Pc 


y 3^2’ 

(2.2) 

Pm 

3772 ’ 

(2.3) 

N = Ina, 

(2.4) 


where a it is the scale factor. It is not hard to see, that for interacting models the autonomous 
system does read as 


= 3x{l + y) - 


Q + 3Hpc{l -|- LOc) 

3H^ 


(2.5) 


and 


y' = 3y{l + y) + 


Pc 

3JL3’ 


( 2 . 6 ) 


where / it is the derivative with respect to N, dot it is the derivative with respect to cosmic 
time, Q it is interaction term. Explicit forms of Eq. (2.5) and Eq. (2.6) are obtained when 
the form of Q is known. To have physically reasonable solutions we should have the following 
constraints 0 < x < 1 and 0 < z < 1. Therefore, we need to find the range of the model 
parameters such, that to have a physically reasonable stable critical points i.e. 0 < Xc < 1 
and 0 < Zc < 1. From the other hand a stable critical point it is an attractor, which we 
are looking for. At the same time we should remember that x and z according to Eq. (1.4), 
Eq. (2.1) and Eq. (2.3) should satisfy to the following constraint 


X + z = 1. 
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(2.7) 







In terms of x and y EoS parameter of the Chaplygin gas does read as 



EoS parameter of the effective flnid does read as 




Pc 

Pc Pm 


( 2 . 8 ) 


(2.9) 


becanse dark matter is considered as a pressnreless flnid, while it is not hard to show that 
the deceleration parameter q does read as 


q = -l-^ = lil + 3y). ( 2 . 10 ) 

In this work we will consider varions cosmological models involving a non linear interaction 
into the darkness of the nniverse. Particnlaiiy, we are interested to consider an interaction Q 
of the following form 

Q = 3Hb(^p+^^, ( 2 . 11 ) 

where p conld be either the energy density of the effective flnid, or the energy density one of 
the components of the effective flnid. p it is prodnct of the energy densities of the components 
of the effective flnid i.e. three possibilities conld be and pcPm- From the next section 

we will start onr analysis. 

3 Phase space analysis 

In this section we will hnd late time attractors of the held eqnations for varions interacting 
varying Chaplygin gas cosmological scenarios. For that we need to determine interaction 
term Q to have Eq. (2.5) and Eq. (2.6). In this work we will pay attention to some hxed 
sign non linear interactions. Non linear interactions having similar strnctnre as here, were 
considered in Ref. [63]. In Ref. [63] varions cosmological models with Chaplygin gas Eq. (1.1) 
have been stndied, inclnding models with some non linear sign changeable interactions. To 
simplify stndy, we will combine an experience known from Physical Literatnre to impose some 
constraints on the model parameters. Namely, for varying Chaplygin gas we will impose that 


0<a<l, (3.1) 

0 < n < 5, (3.2) 

A > 0 and B > 0. (3.3) 

For parameter b describing interaction Q we take 

0 < 6 < 1, (3.4) 


which is the best range snpported from observational data. 6 = 0 does correspond to non 
interacting models. From the other hand, constraints on the deceleration parameter q 



and on the EoS parameter of varying Chaplygin gas 

- 2 < Wc < 0, (3.6) 

will give ns an option to hnd late time attractors for the large scale nniverse with an accelerated 
expansion. It is not exclnded that the phantom divide abont z ~ 0.2 redshift took place in onr 
nniverse, therefore constraint Eq. (3.6) is considered by ns [64] - [66]. Conditional attractors 
are interesting possibilities, becanse depends on the constraints on the model parameters, 
different late time states of the nniverse will be obtained. The goal of this section is to hnd 
and discns stability of the critical points for appropriate cosmological scenarios for 6 different 
cosmological models. 


3.1 Varying Chaplygin gas = Apc - 

We did varions assnmptions allowing ns to stndy cosmological models for the large scale 
nniverse. We start onr analysis of two different cosmological models involving different forms 
of non linear interactions Eq. (2.11) for each of them. We will start from the models where 
varying Chaplygin gas Eq. (1.2) interacts with dark matter. 


3.1.1 Interaction Q = ^Hh (^pc + ^ ^ ^ 

The hrst cosmological model that we wonld like to stndy is a model where interaction between 
varying Chaplygin gas Eq. (1.2) and dark matter is given as 

Q = ^Hh{p, + ^^\ (3.7) 

V Pc P Pm, J 

Snch consideration allowed ns to hnd all real critical points analytically. Obtained solntions 
are presented in Table 1. In this case 3 different physically reasonable critical points are exist 
and only two of them are stable. However only one of them does correspond to the state of 
the large scale nniverse. Stable critical points are E.l.l and Ei.1.3. Critical point E.1.1 is 
late time scaling attractor, becanse 


_ _ 1 - 6 - Vl + (2 - 36)6 

He 2(6-1) ’ ^ 

tends to a constant. For this nniverse the decelerated parameter is q = —1 and Wg// = — 1, 
while the EoS parameter of the varying Chaplygin gas Eq. (1.2) is 


Wc = - 


26 

6- 1 + v'l + (2 - 36)6 


(3.9) 


It is not hard to see that in this case varying Chaplygin gas Eq. (1.2) is a phantom dark 
energy. Disenssed behavior is obtained when the parameters satisfy to Eq.-s (3.1) - (3.3) and 
0 < 6 < 2/3. Critical point Ei.1.3 to be stable when 0 <q;< 1, 0<n<5, 6 = 0 and 
A > 1 + a + n/2 shonld be. In this case a matter will dominate in onr nniverse, since x = 0. 
In snch nniverse accelerated expansion is not possible. Therefore, we left only E.IA as a 
physically reasonable solntion for the large scale nniverse. 
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S.P. X 

y 

fe—l+\/l+(2—36)6 

26 

-1 

pn 1 o A—b+^{A—b)(A+3b) 
2(A-b) 

A(A-b+^J(A-b){A+3b)) 

2{A-b) 

^ A-b-^{A-b)(A+3b) 

2{A-b) 

A{A-b-yJ{A-b){A+3b)) 

2[A-b) 


Table 1. Critical points corresponding to interacting varying Chaplygin gas Eq. (1.2) for interaction 
Q determined via Eq. (3.7). 


S.P. 

X 

y 

E.2.1 

E.2.2 

l+26-v/l+462 

26 

A+2b 

A+b 

-1 

A{a+2b) 

A+b 


Table 2. Critical points corresponding to interacting varying Chaplygin gas Eq. (1.2) for interaction 
Q determined via Eq. (3.10). 


3.1.2 Interaction Q = 3Hb ^pc + 

In the second model the interaction between varying Chaplygin gas Eq. (1.2) and dark matter 
is taken to be 

Q = 3Hb (pc+ ^ . (3.10) 

V Pc Pm J 

In this case physically reasonable two critical points do exist, which are presented in Table 2. 
E.2.1 is a physically reasonable solution, moreover it is late time scaling attractor when 
0<n<5,0<Q;<l,A>0 and 0 < 6 < 2/3. This solution represents a state of the 
universe where varying Chaplygin gas Eq. (1.2) is a phantom dark energy with 


1 + 26 + Vl + 462 

0Jc = - 

2 

This model is free from the cosmological coincidence problem due to 


(3.11) 


Q.m 26 - 1 + Vl + 462 

’■=sf =- Y -■ 

Critical point E.2.2 is physically reasonable solution when 6 = 0, this does mean that x = 1 
and y = A. Therefore, this is the state of the universe where an accelerated expansion is not 
possible, since g > 0. This is the state of the universe where varying Chaplygin gas Eq. (1.2) 
completely dominates in the dynamics of the universe as a usual fluid. As was expected E.2.2 
is not a stable critical point. Fig. (1) does represent phase space portraits for £'.1.2 and £.2.1 
critical points. For a symmetry in plots, we considered — 1 < x < 1 interval, however, we 
should remember that physically reasonable range for x it is (0,1]. 

3.1.3 Interaction Q = 3Hb (^pc + Pm + p^+p 

The other possibility that we are going to consider in this subsection, it is a cosmological 
model where interaction term Q does read as 

Q = 3Hb f Pc -h Pm H--^ • (3.13) 

V Pc + PmJ 
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Figure 1. Phase space portraits for the interacting varying Chaplygin gas Eq. (1.2) models. Left plot 
does represent the model where interaction Q is determined via Eq. (3.7). Right plot does represent 
the model where interaction Q is given via Eq. (3.10). 


S.P. 

X 

y 

E.3A 

2fe-l-|-Vl-462 

2b 

-1 

E.3.2 

A-2b+^/A'2+AAb-4b'^ 

A(A-2b-|-VA2+4Ab-462) 

2{A-b) 

2{A-b) 

E.3.3 

A-2fe-VA2+4Ab-4fe2 

A{A-2b-VA‘^+AAb-Ab^) 

2(A-b) 

2{A-b) 


Table 3. Critical points corresponding to interacting varying Chaplygin gas Eq. (1.2) for Q deter¬ 
mined via Eq. (3.13). 


In this case 3 critical points are physically acceptable (among 4 critical points). They are 
collected into Table. 3. Critical point E.3.1 is a late time scaling solution since 


Qm _ Vl - 462 1 

~ 2(1 - 26) ~ 2’ 


(3.14) 


and does represent the state of the universe, where varying Chaplygin gas Eq. 
phantom dark energy with 

2b 

~ ~26 - 1 + \/l -462' 


(1.2) is a 
(3.15) 


This is possible late time attractor when 0<n<5, 0<6<2/5, 0<a<l and ^4 > 0. 
Critical point E.3.2 is physically reasonable when 6 = 0, and it is unstable, therefore, we will 
not discuss cosmological consequences related to this solution. From the other hand we found 
that E.3.3 is a stable critical point, when 6 = 0 and A > 1 + a + n/2 for 0 < n < 5 and 
0 < a < 1. It is not hard to see that in this case our universe will be in a matter dominating 
state. Accelerated expansion in this case is not possible. 


3.2 Varying Chaplygin gas Pc = Apc - 

re 

In subsections 3.1.1 - 3.1.3 we considered cosmological models where varying Chaplygin gas 
was given via Eq. (1.2). In next three subsections we will consider cosmological models, where 
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S.P. 

EA.l 

EA.2 

EAA 

EAA 


X 


A — b-\-^ (A—fo)(A+36) 
2 {A-b) 

A-b-yJ{A-b){A+1,b) 


2{A-b) 

n+3(6—l)(l+«)—f 
6b(l+a) 

n+3(6—l)(l-j-«)+'f 


6b(l+a) 


y 

A{A-b+^{A-b){A+3b)) 

2{A-b) 

A{A-b-^{A-b){A+3b)) 

2{A-b) 


_1 _|_ n _ 

^ 3(l+a) 


-1 + 


3(l+a) 


Table 4. Critical points corresponding to interacting varying Chaplygin gas Eq. (1.3) for interaction 
Q determined via Eq. (3.7). f = y^(n — 3(1 + q;)(1 + 36))(n — 3(1 + a)(l — 6)). 


another model of varying Chaplygin gas is associated to a part of energy of our universe. 
Namely, we will consider models, where varying Chaplygin gas is given via Eq. (1.3) 


3.2.1 Interaction Q = 3Hb ^pc + ^ ^ ^ 

As in subsection 3.1.1, first of all, we will consider interaction given via Eq. (3.7). In this 
case 4 critical points are obtained (Table 4). To simplify future analysis we can take some 
particular values for the parameters of the model. We will consider a particular example with 
0 < 6 < 1/10, when n = 1. In case of n = 1, when A > b, only critical point EAA is a late 
time scaling attractor with 


r = 


nr 




2 + 3a + 96(1 + a) 


y/(2 + 3a - 36(1 + a))(2 + 3a + 96(1 + a)) j 
Obtained constraints on the parameters provide the deceleration parameter q 

1 

q = -l + 


2{l + a)' 


(3.16) 


(3.17) 


to be 0 < g < — 1. In this case varying Chaplygin gas Eq. (1.3) is a quintessence dark energy 


26(2 + 3a) 

36(1 + a) - 2 - 3a + ^y(2 - 36 - 3a(6 - 1))(2 + 3a + 96(1 + a))' 


(3.18) 


Our study of the general case showed that EAA (only one stable critical point for this model) 
it is a conditional late time scaling attractor. This does mean, that for some constraints on 
the model parameters we can obtain the state of the universe, where varying Chaplygin gas 
Eq. (1.3) will be a quintessence dark energy. From the other hand, EAA could describe the 
state of the universe where varying Chaplygin gas Eq. (1.3) will be a phantom dark energy. 
This is very interesting aspect missing in previously considered models. 


3.2.2 Interaction Q = 3Hb (^pc + 

A cosmological model where varying Chaplygin gas Eq. (1.3) interacts with dark matter 
via interaction Q Eq. (3.10) has three critical points and one of them is a late time scaling 
attractor. Critical points for this model are collected in Table 5. Late time scaling solution 
E.5.2 for some values of parameters could represent a state of the universe, where varying 
Chaplygin gas Eq. (1.3) is a quintessence dark energy (Table 6), while for some other values 
of the parameters we will have a state of the universe, where varying Chaplygin gas Eq. (1.3) 
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S.P. X y 

17 c; 1 A+2b A(A+2b) 

A+b A+b 

17 c: o 3(1+26)(1 +q)— n—r -j . n 

6Wl+a) ^ 3(l+a) 

17 c: o 3(l+2fo)(l+Q)—n+r . 7^ 

6b(l+a) 3(l+a) 


Table 5. Critical points corresponding to interacting varying Chaplygin gas Eq. (1.3) for interaction 
Q determined via Eq. (3.10). f = y^36P^T^r”ap~+~(ri^^'3(rT~ay)^. 


n 

a 

b 

0 < n < 1 

0 < a < 1 

0 < b < 3c,_,_„_,_3 

1 < n < 2 

0 < a < |(2n — 3) 

0<b<^-^^ 

§ < n < 2 

|(2n — 3) < a < 1 

0 < 6 < , 

— 3a+7i+3 

2 < n < 3 

^ <a< i(2n-3) 

0 < 6 < 6«-2ti+6 
^ ^ — 9a+9 

2 < n < 3 

g(2n — 3) < a < 1 

0 < 6 < , 

— 3a+7i+3 

3 < n < 4 

^ < a < 1 

Q <b < 6a-2n+6 
^ ^ — 9a+9 


Table 6. Constraints on the model parameters for the interacting varying Chaplygin gas Eq. (1.3) 
with interaction term Q described. Eq. (3.10). E.5.2 does describe the state of the universe where 
varying Chaplygin gas Eq. (1.3) is a quintessence dark energy. 


n a 

n = 0 0<a<l 

0<n<| 0<a<l 

I < n < 3 5(2n — 3) < a < 1 


0 < 6 < 


3(1+71+3 

n 

3(1+71+3 


< b 

< 


9(1+9 
^ 6a—271+6 
— 9a+9 

r ^ 6a—2n+6 
" — 9a+9 


Table 7. Constraints on the model parameters for the interacting varying Chaplygin gas Eq. (1.3) 
with interaction term Q described. Eq. (3.10). E.5.2 does describe the state of the universe where 
varying Chaplygin gas Eq. (1.3) is a phantom dark energy. 


is a phantom dark energy (Table 7). To finish this subsection we provide expressions for the 
deceleration parameter q, Uc and r 


and 


where r 


q = -l + 


n 

2(1 +a)’ 


Wc 


n 

^e//- -1+3(1^^), 

3(1 + 26)(1 + a) — n + f 

6(1 + a) 


3(1 — 26)(1 + a) — n — r 
^ “ 2(n-3(l + a)) ’ 

Y^3662(1 + q;) 2 + (n — 3(1 + a))^, for E.5.2 attractor. 


(3.19) 

(3.20) 

(3.21) 

(3.22) 


3.2.3 Interaction Q = 3Hb (^pc + Pm + ^ ^ ^ 

The last model that we will consider in this work it is a model where the interaction between 
varying Chaplygin gas Eq. (1.3) and dark matter is given via Eq. (3.13). For this model 4 real 
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S.P. 

X 

y 

E.6.1 

A-2fe+A/yl2+4Ab-462 

A(A-2b+VA2+4Ab-462) 

2{A-b) 

2{A-b) 

E.6.2 

A-2fe-VA2+4Ab-4fe2 

A( A-2b+VA2 +4Ab-462) 

2(A-b) 

2{A-b) 

E.6.3 

n-|-3(—l-i-2b) (l+o)—r 

1 1 ” 

6fe(l-|-o) 

^ ' 3(l+a) 

E.6A 

n-|-3( —l-|-2b)(l-|-Q:)-|-r 

1 1 

6fe(l-|-o) 

^ ' 3(l+a) 


Table 8. Critical points corresponding to interacting varying Chaplygin gas Eq. (1.3) for interaction 
Q determined via Eq. (3.13). f = + a)^ + (n — 3(1 + a))^. 


n 

a 

b 

n = 0 

0 < n < 1 

1 < n < 3 

0 < a < 1 

0 < a < 1 
^(2n — 3) < a < 1 

0 ^ h 18a2+36Q.+i8 

^ ^ 45q:2+90q+45 

—n2+3Q:n+3n ^ u ^ 6a—2n+6 

9a2+18Q-|-n2+9 ^ — 15a+15 

—n2+3on+3n ^ u ^ 6a—2n+6 

9Q;^+18QH-n^H-9 — 15 q:+15 


Table 9. Constraints on the model parameters for the interacting varying Chaplygin gas Eq. (1.3) 
with interaction term Q described. Eq. (3.13). E.6A does describe the state of the universe where 
varying Chaplygin gas Eq. (1.3) is a phantom dark energy. 


critical points do exist, which are presented in Table 8. However, only 
scaling attractor. Moreover, this attractor it is a conditional attractor, 
the following expressions 


q = -l + 


n 

2(1 +a)’ 


E.QA is a late time 
It is easy to obtain 

(3.23) 




3(l + a)’ 


(3.24) 


and 


26(ra-3(l + a)) 
n + 3(—1 + 26)(1 + a) + f 


1 f 

2 “ 2(n-3(-l + 26)(l + a))’ 


(3.25) 


(3.26) 


where r = -^—365^(1 + a)^ + (n — 3(1 + a))^, for E.GA solntion. In Table 9 constraints on 
the model parameters are collected, when attractor E.6A does describe a nniverse where 
varying Chaplygin gas Eq. (1.3) is a phantom dark energy. It is also possible, that for some 
valnes of the model parameters varying Chaplygin gas Eq. (1.3) is a qnintessence dark energy. 
It is mathematically very easy problem and we omitted appropriate discnssion on this issne 
to save a place. 


4 Discussion 

In this work we have considered 6 different cosmological scenarios. Three different strnctnres 
for the interaction term Q are nsed. Chaplyging gas (and varions modifications of it) is still 
very active snbject for stndy in cosmology, becanse it is a joint model of dark energy and 
dark matter. One of the modifications of it, it is varying Chaplygin gas where a phenomeno¬ 
logical assnmption abont B (or A) parameter is nsed. We considered two models where two 
parametrizations of B are taken into acconnt. Assnming the strnctnre of the darkness of the 
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universe and the form of the interaction term Q we performed phase space analysis. For all 
cases late time scaling attractors are found and appropriate constraints on the model parame¬ 
ters are obtained. For the first type cosmological models late time scaling attractors describe 
the state of the universe where varying Chaplygin gas it is a phantom dark energy. In the 
second models attractors are conditional late time scaling attractors. This does mean that for 
some values of the model parameters varying Chaplygin gas will be a phantom dark energy, 
while for some cases it will be a quintessence dark energy. An appropriate discrimination of 
discussed issues are presented in this article to give a right motivation for other works. We 
hope that very soon we will report about our new results on this topic. 
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